A  MATHEMATICAL  MODEL  FOR  THE  STUDY  OF  HEMORRHAGIC 
SHOCK  AND  FLUID  RESUSCITATION: 
TRANSCAPILLARY  EXCHANGE 


Tammy  J.  Doherty 


Division  of  Military  Trauma  Research 


LETTERMAN  ARMY  INSTITUTE  OE  RESEARCH  PRESIDIO  OF  SAN  FRANCISCO  CALIFORNIA  94129 


A  MATHEMATICAL  MODEL  FOR  THE  STUDY  OF  HEMORRHAGIC 
SHOCK  AND  FLUID  RESUSCITATION:  TRANSCAPILLARY  EXCHANGE 
-  Tammy  J.  Doherty 


This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  to  the 
originator. 

Citation  of  trade  names  in  this  report  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such  items. 


This  material  has  been 
reviewed  by  Letterman  Army 
Institute  of  Research  and  there 
is  no  objection  to  its 
presentation  and/or  publication. 
The  opinions  or  assertions 
contained  herein  are  the  private 
views  of  the  author(s)  and  are 
not  to  be  construed  as  official 
nor  as  reflecting  the  views  of 
the  Department  of  the  Army  or 
the  Department  of  Defense. 
(AR  360-5) 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMB  No  0704-01 81 


i*  REPORT  security  Classification 

Unclassified 

la  SECURITY  CLASSIFICATION  authority 

2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 

Institute  Report  No.  477 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

6b.  OFFICE  SYMBOL 

Letterman  Army  Institute 

(If  applicable) 

of  Research 

SGRD-ULT-M 

6c.  ADDRESS  (City,  State,  and  HP  Code) 

Letterman  Army  Institute  of  Research 

Division  of  Military  Trauma  Research 

Presidio  of  San  Francisco,  CA  94129-6800 

8a.  NAME  OF  FUNDING /SPONSORING 

8b  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

Div.  Military  Trauma  Res. 

SGRD-ULT-M 

8c  ADDRESS  (City,  State,  and  ZIP  Code) 

lb  RESTRICTIVE  markings 


3  DiSTRiBuT.ON/ availability  Of  REPORT 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a  NAME  Of  MONITORING  ORGANIZATION 
U.S.  Army  Medical  Research  and  Development  Command 


7D  ADDRESS  (Ciry.  Sure.  and  ztf  Code) 
Fort  Detrick 

Fredrick,  MD  21701-5012 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO  NO 

61101A 


11.  TITLE  (include  Security  Classification) 

(U)  A  Mathematical  Model  for  the  Study  of  Hemorrhagic  Shock  and  Fluid  Resuscitation:  Transcapillary  Exchange 


12.  PERSONAL  AUTHOR(S) 
Tammv  J.  Doher 


13*  TYPE  OF  REPORT 

Final 


T  3b.  TIME  COVERED 
FROM  1QRQ  TO  Itm 


14.  DATE  OF  REPORT  (Year,  Month.  Day)  |1S.  PAGE  COUNT 
March  1993  I  32 


17.  COSATI  COOES  |  18.  SUBJECT  TERMS  (Continue  on  reverie  if  necessary  and  identify  by  block  number) 

FIELD  I  GROUP  I  SUB-GROUP  I  ,  .  ,  »  „  ...  _ 


model,  hemorrhage,  fluid  resuscitation,  transcapillary  exchange 


19  ABSTRACT  (Continue  on  rev ene  if  neceisery  end  identify  by  block  number)  - \ 

Mathematical  models  used  to  predict  plasma  volume  responses  to  hemorrhage  and  fluid  resuscitation  must  consider  both 
fluid  and  solute  exchange  across  the  capillary  wall.  In  this  paper,  we  review  four  models  of  transcapillary  exchange  that  might 
be  incorporated  into  a  comparlmental  model  for  plasma  volume  prediction.  The  simplest  of  the  four  transcapillary  flux  models, 
attributed  to  Kedem  and  Katchalsky  (1958),  assumes  a  discrete  capillary  wall  structure,  ideal  solutions,  and  capillary  wall  j 
homogeneity.  The  remaining  models  become  progressively  more  complex  by  successive  elimination  of  these  three  assumptions,  j 
To  assess  their  effects  on  whole-body  plasma  volume  predictions,  the  four  models  of  transcapillary  exchange  were  incorporated 
into  a  simple  two-compartment  (blood  and  interstitial  space)  model.  Next,  plasma  volume  predictions  were  generated  for  the 
following  simulation  conditions:  1)  instantaneous  25%  loss  of  blood  volume,  2)  instantaneous  50%  increase  in  plasma  NaCl 
concentration,  3)  instantaneous  50%  increase  in  plasma  albumin  concentration,  and  4)  instantaneous  25%  loss  of  blood  volume 
with  a  50%  increase  in  plasma  NaCl  concentration.  Plasma  volume  predictions,  generated  by  the  four  models,  were 
indistinguishable  over  these  simulation  conditions.  From  these  results,  we  conclude  that  the  simplest  of  the  transcapillary  flux 
models  may  be  used  for  predicting  plasma  volume  responses  to  hemorrhage  and  resuscitative  fluid  administration. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
03 UNCLASSIFIEDrtJNUMlTED  □  SAME  AS  RPT.  □  OTIC  USERS 


21  ABSTRACT  SECURITY  CLASSIFICATION 


22b  TELEPHONE  (include  Area  Code )  1  22c  OFFICE  SYMBOL 


iWTA 


DD  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


ABSTRACT 


Mathematical  models  used  to  predict  plasma  volume  responses  to 
hemorrhage  and  fluid  resuscitation  must  consider  both  fluid  and  solute 
exchange  across  the  capillary  wall.  In  this  paper,  we  review  four  models  of 
transcapillary  exchange  that  might  be  incorporated  into  a  compartmental 
model  for  plasma  volume  prediction.  The  simplest  of  the  four 
transcapillary  flux  models,  attributed  to  Kedem  and  Katchalsky  (1958), 
assumes  a  discrete  capillary  wall  structure,  ideal  solutions,  and  capillary 
wall  homogeneity.  The  remaining  models  become  progressively  more 
complex  by  successive  elimination  of  these  three  assumptions.  To  assess 
their  effects  on  whole-body  plasma  volume  predictions,  the  four  models  of 
transcapillary  exchange  were  incorporated  into  a  simple  two-compartment 
(blood  and  interstitial  space)  model.  Next,  plasma  volume  predictions 
were  generated  for  the  following  simulation  conditions:  1)  instantaneous 
25%  loss  cf  blood  volume,  2)  instantaneous  50%  increase  in  plasma  NaCl 
concentration,  3)  instantaneous  50%  increase  in  plasma  albumin 
concentration,  and  4)  instantaneous  25%  loss  of  blood  volume  with  a  50% 
increase  in  plasma  NaCl  concentration.  Plasma  volume  predictions, 
generated  by  the  four  models,  were  indistinguishable  over  these  simulation 
conditions.  From  these  results,  we  conclude  that  the  simplest  of  the 
transcapillary  flux  models  may  be  used  for  predicting  plasma  volume 
responses  to  hemorrhage  and  resuscitative  fluid  administration. 

keywords:  model,  hemorrhage,  fluid  resuscitation,  transcapillary  exchange 


......  3 


DOHERTY  -  1 


A  Mathematical  Model  for  the  Study  of  Hemorrhagic  Shock  and  Fluid 
Resuscitation:  Transcapillary  Exchange  --  Tammy  J.  Doherty 


INTRODUCTION 

One  of  the  most  important  recovery  mechanisms  following  loss  of 
blood  volume  is  the  transfer  of  fluid  across  the  capillary  wall  from  the 
interstitium  (subscript  I)  to  the  capillary  (subscript  C).  According  to 
Starling1,  the  forces  that  govern  transcapillary  fluid  movement  are  the 
hydrostatic  (P)  and  osmotic  (II)  pressure  differences  across  the  capillary 
wall.  Landis2  showed  that  the  rate  of  fluid  movement  from  capillary  to 
interstitium  (Jv)  is  approximately  proportional  to  the  difference  between 
these  capillary  (subscript  C)  and  interstitial  (subscript  I)  pressure 
differences  (Box  1): 

Jy  =  k[(Pc  -  p7)  -  (nc  -  n,)]  (1) 

When  capillary  hydrostatic  pressure  falls  due  to  loss  of  blood  volume,  fluid 
movement  into  the  capillary  space  increases,  helping  restore  lost  volume. 
This  "autoinfusion"  dilutes  the  blood  (lowering  osmotic  pressure)  and 
causes  capillary  hydrostatic  pressure  to  increase.  Eventually,  a  new  balance 
of  hydrostatic  and  osmotic  forces  is  reached  and  transcapillary  flux  returns 
to  normal  levels. 

Pre-hospital  management  of  hemorrhagic  shock  consists  of:  1) 
establishing  hemostasis,  and  2)  replacing  lost  blood  with  a  fluid  substitute. 
Because  whole  blood  is  oftentimes  unavailable,  lactated  Ringer’s  solution, 
which  has  the  same  osmolarity  as  normal  plasma,  is  frequently  used  as  a 
fluid  substitute.  Fluid  and  solute  components  of  lactated  Ringers  solution 
distribute  throughout  the  interstitial  as  well  as  the  vascular  spaces. 
Consequently,  it  takes  three  or  more  liters  of  lactated  Ringers  solution  to 
replace  one  liter  of  shed  blood.  When  large  quantities  of  blood  or 
isosmotic  fluid  are  unavailable,  a  small  volume  of  a  hyperosmotic  solution 
may  be  used  to  expand  plasma  volume  using  the  body’s  own  autoinfusion 
mechanism.  Capillary  permeability  to  the  solutes  in  the  resuscitative  fluid 
determines  the  duration  of  the  osmotic  pressure  differential  across  the 
capillary  wall.  A  hyperosmotic  solution  made  up  of  highly  permeable 
solutes  (e.g.,  NaCl)  will  have  only  a  short  term  effect  compared  to  solutions 
made  up  of  less  permeable  solutes  (e.g.,  albumin,  dextran,  hetastarch), 
which  stay  in  the  blood  space  longer.  For  these  reasons,  the  dynamics  of 
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both  solute  and  solvent  must  be  considered  when  interpreting 
cardiovascular  responses  to  hemorrhage  and  fluid  resuscitation. 

The  Landis-Starling  Equation  (Equation  1)  does  not  account  for  the 
permeability  of  solutes.  Transport  equations  that  account  for  permeability 
of  both  solute  and  solvent  generally  fall  into  one  of  two  classes: 
phenomenologic  equations  based  on  irreversible  thermodynamics,  or 
mechanical  (hydrodyn*:  mic)  equations  based  on  either  "pore"  or  "matrix" 
theory.  The  models  in  these  two  classes  differ  mainly  with  respect  to  the 
definition  of  model  coefficients:  the  forms  of  the  transport  equations  are 
otherwise  similar.  Consequently,  models  may  be  characterized  by  the 
simplifying  assumptions  used  in  model  development,  rather  than  in  the 
theory  used  to  develop  equations. 

It  is  difficult  to  determine,  on  the  basis  of  model  descriptions,  which 
model  or  set  of  model  assumptions  is  necessary  or  sufficient  for  application 
to  the  problem  of  hemorrhage  and  fluid  resuscitation.  The  simplest 
transport  model,  attributed  to  Kedem  and  Katchalsky3  has  been  criticized 
because  it  incorporates  assumptions  of  a  thin,  homogeneous  capillary  wall, 
and  ideal  solutions  made  up  of  small  particles.  In  this  paper,  we  review 
the  Kedem-Katchalsky  model  as  well  as  other,  less  restrictive  models.  The 
purpose  of  this  endeavor  is  to  identify  the  simplest  model  for  predicting 
whole-body,  net  transcapillary  fluid  and  solute  exchange  in  response  to 
hemorrhage  and/or  resuscitation  fluid  administration. 


MODEL  DESCRIPTIONS 


Phenomenologic  models  based  on  irreversible  thermodynamics  were 
first  introduced  by  Kedem  and  Katchalsky3  in  1958.  The  model  they 
developed  assumes  a  homogeneous,  discrete  membrane  structure,  and  ideal 
solutions  composed  of  a  single  type  of  small  solute  particle.  Several 
models  have  been  developed  to  address  perceived  shortcomings  in  the 
Kedem-Katchalsky  model.  These  include  a  model  developed  by  Patlak4  in 
1963,  which  uses  a  continuous  membrane  model,  as  well  as  Katz’  model5, 
developed  in  1985,  which  assumes  non-ideal  solutions.  The  basic  pore 
theory  model,  described  by  Curry6,  also  incorporates  a  continuous 
membrane  model  and  non-ideal  solutions.  In  addition,  the  pore  theory 
model  readily  accounts  for  large  solute  molecules  and  multiple  pathways 
for  solute  and  solvent  transport.  To  better  understand  the  effects  of  model 
assumptions  on  flux  equations,  three  models  based  on  irreversible 
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thermodynamics  (i.e.,  the  Kedem-Katchalsky,  Patlak,  and  Katz  models),  as 
well  as  the  basic  pore  theory  model,  are  reviewed  below. 


The  Kedem-Katchalsky  Equations 

In  models  based  on  irreversible  thermodynamics,  the  rate  of  entropy 
production  (dS/dt)  on  one  side  of  the  membrane  is  defined  in  terms  of 
chemical  potentials  (y).  This  equation  is  then  used  to  generate  coupled 
total  volumetric  and  solute  flux  equations.  Models  differ  with  respect  to 
the  form  of  the  entropy  production  equation  (discrete  versus  continuous) 
and  with  respect  to  simplifying  assumptions  used  in  estimating  chemical 
potentials.  In  the  Kedem-Katchalsky  model,  'he  rate  of  entropy  production 
for  the  interstitial  compartment  is  given  by: 


dSj 

~dt 


— I  =  —Au 

t.  ry,  *  y 


dN 


dt 


Ms- 


dN. 


S,1 


dt 


(2) 


where  T  is  temperature,  t  is  time,  N  is  the  number  of  moles  of  solute 
(subscript  s)  or  water  (subscript  w),  and  the  A  operator  represents  capillary- 
interstitial  differences.  For  convenience,  the  variable  $  is  defined  as: 

j*  dSj 

= - =Au/t+Aun  w) 

A  dt  “w  w  r s  s 


where  A  is  the  membrane  surface  area  and  fy  is  defined  by: 


A 


(4) 


For  an  ideal  solution,  the  chemical  potential  difference  (ti)  is  defined  in 
terms  of  hydrostatic  pressures  (P),  molar  volumes  (V),  and  partial  molar 
fractions  (7): 

A  \i  =  VAP  +  RT  Aln(y)  (5) 


where  R  is  the  universal  gas  constant  and  the  partial  molar  solute  fraction 
is  defined  by: 


Y-  = 


nS  +  nw 


(6) 


Assuming  a  dilute,  binary  solution,  (i.e.,  n^r^  and  iij  +  ivm,,),  7S  may  be 
approximated  by  the  term  (n^/nj  and  solute  concentration  (Cs)  by  the 
term  (nj/VjT,,).  Using  these  approximations,  Aln(Cs)  may  be  substituted 
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for  Aln(ts)  in  Equation  5,  and  the  solute  chemical  potential  difference 
across  the  membrane  becomes: 

Ap,  =  V^AP  +  RT  AlnfC,)  (7) 

In  their  paper3 ,  Kedem  and  Katchalsky  rewrite  Equation  7  in  the  following 
form: 


A  \is  =  VSA  P  +  RT 


(8) 


where  C  is  defined  by: 


Aln(C,) 


(9) 


To  determine  the  corresponding  equation  for  the  water  or  solvent  chemical 
potential  difference,  ln(7w)  is  represented  by  ln(l-7s)  which,  for  small  7S,  is 
approximately  equal  to  -7S.  Again  assuming  dilute,  binary  solutions  (i.e., 
Hs+Hw-Hw  and  Q^/V^),  the  expression  for  the  water  chemical  potential 
difference  becomes: 

A ]iw  =  VWA P  -  RTVWACS  (10) 


Substituting  the  definitions  of  A/xs  and  from  Equations  8  and  10  into 
Equation  3  and  rearranging  terms  yields: 


ns  .  „ 

-r  -  KK 

\C* 


RTAC ' 


(11) 


Equation  11  may  be  simplified  by  adding  and  subtracting  the  term 
r^VjRTAC,,  by  substituting  Jv  for  the  term  (r^V.+r^VJ,  and  substituting 
All  for  the  term  (RTACS): 


#  =  JV(AP  -  AH)  ♦  JsV\  \  + 


1  ^ 


csK ) 


AH 


(12) 


By  the  rules  of  irreversible  thermodynamics7,  any  equation  for  the  rate  of 
entropy  generation  of  the  form  of  Equation  12  (i.e.,  sums  of  fluxes 
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multiplied  by  their  own  driving  forces)  may  be  recast  in  the  form  of  exact 
coupled  flux  equations  as  follows: 


Jv  =  Ln(AP  -  AD)  ♦  Ln 


1  + 


AH 


(13) 


JsVs  -  ^(A P  -  AD)  ♦  Lj\* 


C  V 

v  s  ) 


ah 


(14) 


An  expression  for  (AP-AII)  may  be  obtained  from  Equation  13  and 
substituted  into  Equation  14  to  yield: 


JsK  - 


'21 


'11  1 


4,  ^12 


1  + 


C*K, 


AD 


+  L  22  i  + 


CsKj 


AH 


(15) 


By  assuming  that  (Onsager’s  reciprocity  law),  that  AII  =  RTACS 

(van’t  Hoffs  law  for  ideal  solutions),  and  that  CsVs»0  for  small  particles  in 
dilute  solution,  Equations  13  and  15,  after  rearranging  terms,  become: 

Jv  =  *V[AP  -  oAH]  (16) 


Js  =  *,A C,  ♦  (1  ~o)JvCs 

where: 

kv  -  Ln 


(17) 

(18) 


o  =  l 


'12 


'11  1 


1  + 


CsK; 


(19) 
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RT 


L72 


Ln, 


1  + 


1  ^ 


C  V 

r  s  ! 


(20) 


Equation  16  is  similar  to  the  Landis-Starling  Equation  (Equation  1)  but 
includes  an  "osmotic  reflection  coefficient"  (0),  defined  as  the  fraction  of 
the  osmotic  pressure,  that  would  be  manifest  across  an  ideal, 
semipermeable  membrane,  that  is  manifest  across  the  non-ideal  capillary 
wall8.  From  this  definition,  0=  1  when  the  membrane  is  impermeable  to 
solute,  and  a«0  when  the  membrane  is  highly  permeable  to  solute.  A 
summary  of  the  Kedem-Katchalsky  Equations,  extended  for  a  system  of 
multiple,  non-interacting  solutes,  is  presented  in  Box  2. 


The  Patlak  (nonlinear)  Equations 

A  major  assumption  used  in  the  Kedem-Katchalsky  formulation  is 
that  the  capillary  wall  is  thin  enough  to  represent  pressure  and 
concentration  gradients  in  discrete  form.  This  assumption  has  been  shown 
to  be  invalid  unless  the  rate  of  solvent  flux  is  very  small4  911.  Patlak4  was 
the  first  to  reformulate  the  solute  and  solvent  flux  equations  using  a 
continuous  membrane  model: 


Jv  =  *vw 


-  0— 
dx  dx 


(21) 


Js  =  ~kt»-r  +  O-oVvC,  {  } 

s  dx 

where  x  is  a  point  along  the  transport  pathway,  and  w  is  the  total  distance 
across  the  capillary/interstitial  membrane.  Assuming  that  k^,  Iq,  and  o  are 
constant  (i.e.,  independent  of  Cs  and  x),  and  that  a  steady  state  is  obtained 
(i.e.,  Js  and  Jv  are  independent  of  time  and  distance,  x),  then  integration  of 
Equations  21Vand  22  leads  to  the  following  solutions  for  Jv  and  Js: 

Jv  -  kv(dF  -  oAII)  (2^) 


Js  -  (i-o K 


Pe) 

(1-e-*) 


(24) 


where  Pe,  a  form  of  the  Peclet  number,  is  defined  as: 
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Pe 


(1  ~o)Jv 

K 


(25) 


Equations  23-25,  extended  to  a  system  of  multiple,  non-interacting  solutes, 
are  presented  in  Box  3. 


Equations  for  non-ideal  solutions 

Both  the  Kedem-Katchalsky  and  Patlak  equations  assume  that  the 
solutions  on  the  two  sides  of  the  membrane  are  ideal.  Katz5  rederived  the 
transport  equations  using  a  continuous  membrane  model,  assuming  non¬ 
ideal  solutions.  As  in  the  Kedem-Katchalsky  derivation,  Katz  begins  his 
derivation  by  defining  a  variable  $  which  is  the  rate  of  entropy  production 
multiplied  by  temperature  and  divided  by  membrane  area.  Because  Katz 
assumes  a  continuous  membrane,  the  expression  for  $  differs  somewhat 
from  the  expression  used  by  Kedem  and  Katchalsky  (see  Equation  3): 


_  TdS  _  .  ^  +  .  d\is 
"  A  dt  "  ”w  dx  +  n*  dx 


(26) 


where  n  is  defined  as  in  Equation  4.  Katz  uses  a  general  expression  for 
chemical  potential: 

p  =  VP  +  RT\n(a)  (27) 


in  which  a  is  the  "chemical  species  activity",  rather  than  the  more  restrictive 
definition  used  by  Kedem  and  Katchalsky  (Equations  5-10).  Taking  the 
differential  of  both  sides  of  Equation  27,  it  is  possible  to  write: 


dp. 

-  K- 

d\n(as) 

♦  RT  * 

(28) 

dx 

s  dx 

dx 

dpw 

m  v  dP 

dln(aj 

+  RT 

(29) 

dx 

w  dx 

dx 

Substituting  these  definitions  of  d ps/dx  and  dpw/dx  into  Equation  26  and 
rearranging  terms  yields: 


dP 

dx 


.  d\n(as)  .  dln(aw)>| 


dx 


+  n. 


dx 


RT 


(30) 
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Adding  and  subtracting  the  term: 

nsVsRT  dln(aw) 
K  dx 


(31) 


to  and  from  the  right  side  of  equation  30,  substituting  Jv  for  the  term 
(r^Vj+r^.VJ,  substituting  Js  for  the  term  (n* ),  and  regrouping  terms  yields: 


RT  d\n{aw) 
+  Vw  dx 


+  J  V 

^sT  s 


RT  dWas) 
Vs  dx 


RT  d\n{aw))  (32) 
dx 


In  an  isothermal,  noncontractile,  isoelectric  system  with  Nw  moles  of  water 
and  Ns  moles  of  solute,  the  Gibbs-Duhem  equation5  states: 

-VdP  +  =  0  (33) 


or,  in  the  present  case: 


-  {nwVw+naV)dP  +  nw[VwdP  +  RTd\n(aw))  *  ns[VsdP  +  RTd]n(as)) 


The  dP  terms  in  Equation  34  drop  out,  leaving  an  expression  for  dln(as)  as 
a  function  of  dln(aw): 


</ln(a4)  = 


-  —  dln(aw) 
n. 


(35) 


Substituting  this  expression  for  dln(as)  into  Equation  32  and  rearranging 
terms  yields: 


> 

RT  dln(awy 

-RT  dln(flw)'| 

dx 

K  dx. 

*  T  ■'X.l 

K  dx 

(36) 


Katz  defines  the  osmotic  pressure  differential  as  the  hydrostatic  pressure 
gradient  that  must  be  applied  to  a  system  to  maintain  that  system  in 
equilibrium  when  a  chemical  potential  gradient  exists5: 
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dU  _  RT  d]n(aw) 

dx  ~  ~  Vw  dx 


Using  this  definition,  we  may  rewrite  the  expression  for  $  (Equation  36)  in 
terms  of  II: 

•  (38) 

"U  dx)  1  •{  * ,V,jdx 


By  the  rules  of  irreversible  thermodynamics7,  any  equation  for  the  rate  of 
entropy  generation  of  the  form  of  Equation  38  (i.e.,  sums  of  fluxes 
multiplied  by  their  own  driving  forces)  may  be  recast  in  the  form  of  exact 
coupled  flux  equations  as  follows: 


’  “1*  dx  j  12  (  naVx  )dl 


*  *  ^ l\dx  dx)  nsVs  dx 


Using  the  same  procedures  as  in  the  Kedem-Katchalsky  derivation,  these 
equations  may  be  simplified  to: 

j.  **[£-„"]  Hi 

•  ’  dx  dx 


J '■  *  *  a>J- 


where: 


*  -in 

% 

w 
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o  =  l- 


'n 


1  + 


n 

IM 

w 

w 

n. 

V. 

s 

*  / 

(44) 


1 

a  =  - 

1 


( 


hi 


«2  = 


r*\ 

'12 


-n; 


+ 

l  *sK) 


(45) 


*02 

hK 


(46) 


The  solution  for  Equation  41  is  straightforward  and  is  identical  to  the 
Patlak  solution  (Equation  23).  To  solve  Equation  42,  Katz5  rewrites  the 
equation  in  terms  of  solute  concentration,  using  a  third  order  polynomial  to 
express  osmotic  pressure  as  a  function  of  solute  concentration: 


n  ~  aiCf  +  &iCs  +  &jCs 


(47) 


Taking  the  derivative  of  both  sides  of  Equation  47  with  respect  to  Cs 
results  in  the  following  relationship  between  dll  and  dCs: 

dU  =  (o,  +2^0,  ♦  3 03  C?)dCs  W) 


Katz  substitutes  this  expression  for  dll  into  Equation  42,  resulting  in  the 
following  expression  for  Js  in  terms  of  Cs: 

dC 

Js  =  «1W(°1  +  2<h  Cs  +  3°3C*)-^£  +  a2JV 

to  solve  this  equation,  Katz  assumes  that  the  o,  term  is  concentration- 
dependent: 


a,  =  1  +  —  (50) 

c* 

and  that  the  a2  term  is  concentration-independent.  Using  these 
assumptions,  it  is  possible  to  solve  Equation  49: 
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Js  =  [60A(ln(Q)  +  A(C,)  +  ^A(C?)  +  fc3A(C43)]  ♦  a2Jv  (51) 

where,  the  b;  terms  are  functions  of  ft3  and  a;  terms,  and  as  before,  the 
symbol  A  represents  capillary-interstitial  differences. 


Unfortunately,  the  assumptions  of  constant  a,  and  concentration- 
dependent  a2,  used  to  solve  Equation  49,  contradict  accepted  notions  of 
concentration  independence  for  the  coefficient  of  diffusion  (related  to  a:) 
and  concentration-dependence  of  the  convective  transport  coefficient  (a2), 
which  is  usually  represented  by  the  product  of  the  solute  concentration 
(Cs),  and  a  "slip"  coefficient  (e.g.,  1  -a  in  the  Kedem-Katchalsky  solution). 
Because  there  is  no  apparent  basis  for  the  assumptions  used  by  Katz,  an 
alternative  solution  is  presented  below. 


To  begin,  an  expression  for  Js  is  derived  from  Equations  39  and  40 
using  the  coefficients  w  and  a: 


where: 


J. 


-  uw 


dn 

dx 


♦  (l-o)  JVC, 


0) 


-Lfe-Jf, 


(52) 


(53) 


and  a  is  defined  as  in  Equation  44.  Substituting  the  expression  for  dll  from 
Equation  48  into  Equation  52,  and  solving  the  resulting  differential 
equation  yields  the  following  equation  relating  solute  concentration,  solute 
flux,  and  total  volumetric  flux: 


k 


2a1Js  3o3y2 
b  +  b2  , 


+  2a2^CJ  +  3a3&Cs 


CO 

(54) 


where: 
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b  =  (l-o)  Jv 


(55) 


—  c  +  c 

c  =  -L*J-LL  (56) 

*  2 

An  approximate  solution  for  Js  from  Equation  54  may  be  obtained  using 
numerical  search  techniques.  The  new  non-ideal  flux  equations  (i.e., 
Equations  41  and  54),  extended  to  a  system  of  multiple,  non-interacting 
solutes,  are  presented  in  Box  4. 


Pore  Theory 

The  first  discussion  of  pore  theory  is  attributed  to  J.R. 
Pappenheimer12.  Since  that  time,  numerous  descriptions  of  hydrodynamic 
flow  through  pores  traversing  the  capillary  wall  have  been  published.  In 
this  review,  we  assume  that  both  water  and  solutes  (approximated  by  rigid 
spheres)  are  transported  across  the  capillary  wall  through  cylindrical  pores. 
The  conceptual  model  by  Anderson  and  Malone13  (Figure  1)  is  used  to 
estimate  both  solvent  and  solute  flux  through  these  pores.  The  distance 
from  the  center  of  the  pore  to  the  wall  is  rp,  and  the  radius  of  the  solute 
molecule  is  rs.  To  determine  the  solute  concentration  at  any  location 
within  the  pore,  the  centers  of  the  solute  molecules,  within  the  region 
bounded  by  x+Ax,  are  counted.  Because  the  solute  molecules  are  of  finite 
size,  the  centers  of  the  molecules  are  prohibited  from  the  region  adjacent 
to  the  pore  wall.  Thus,  the  capillary  pore  may  be  divided  into  a  zone  near 
the  pore  axis  (i.e.,  the  "core  region"),  in  which  solute  is  free  to  move 
(0<r<rp-rs),  and  a  zone  near  the  pore  wall  (the  "exclusion  region"),  in  which 
the  solute  concentration  is  zero  (rp-rs  <  r  <  rp). 

A  key  assumption  in  this  model  is  that  of  radial  equilibrium. 
Hydrostatic  and  osmotic  pressure  gradients  are  assumed  not  to  exist  in  the 
radial  (r)  direction  except  across  the  imaginary  boundary  at  r=rp-rs.  Thus, 
the  hydrostatic  and  osmotic  pressures  at  any  point  in  the  "core"  region,  at 
distance  x  along  the  pore,  may  be  expressed  by  the  hydrostatic  and  osmotic 
pressures  at  the  pore  axis  (P(x)  and  II(x),  respectively).  Because  the  solute 
concentration  in  the  exclusion  region  is  zero,  the  osmotic  pressure  in  that 
region  must  also  be  zero.  Across  the  imaginary  boundary  at  r  =  rp-rs,  the 
driving  force  for  solvent  flow  is  the  difference  between  hydrostatic  and 
osmotic  pressures.  Therefore,  the  hydrostatic  pressure  in  the  exclusion 
region  equals  P(x)-Il(x). 


DOHERTY  -  13 


Using  these  results,  i.e.,  that  the  hydrostatic  pressure  in  the  core 
region  is  P(x),  and  the  hydrostatic  pressure  in  the  exclusion  region  is  [P(x)  - 
n(x)],  it  is  possible  to  write  simplified  Navier-Stokes  Equations  for  the  core 
and  exclusion  regions,  separately: 


where  r?  is  fluid  viscosity.  By  separating  Equation  58  into  hydrostatic 
(subscript  P)  and  osmotic  (subscript  Ii)  pressure  components: 


dP  r\  d  rdvr 
dx  r  dr  dr t 


V 


r,*r*rP 


(59) 


dn  s 

dx  r  dr {  dr)' 


-r<r*rp 


(60) 


it  is  possible  to  generate  equations  (using  Equations  57,  59,  and  60)  that 
describe  the  effects  of  hydrostatic  pressure  over  the  entire  pore  region 
(Equation  61),  and  osmotic  pressure  in  the  exclusion  region  (Equation  62): 


dP  m  rjjrf 
dx  r  dr  \  dr  , 


Os >rsrp 


(61) 


dH  =  JQ. 

dx  r  dry  dr  t 


r-r  <rir„ 

p  s  p 


(62) 


The  solute  velocity  due  to  hydrostatic  pressure  may  be  obtained  by 
integrating  Equation  61  twice  with  respect  to  r,  applying  the  boundary 
conditions:  1)  dvP/dr=0  at  r=0,  and  2)  vP  =  0  at  r  =  rp: 
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v/> 


(^~rl)dP 

4rj  dx 


(63) 


Similarly,  velocity  due  to  osmotic  pressure  in  the  exclusion  region  may  be 
obtained  by  integrating  Equation  62  twice  with  respect  to  r,  applying  the 
boundary  conditions:  1)  vn  =  0  at  r  =  rp-rt  and  2)  vn  =  0  at  r  =  rp: 


vn  “ 


-(r2~rp)dH 
4q  dx 


(.rP~rs)2  dll 
2i\  dx 


l  In  O’) 


ln(^)] 


(64) 


The  solvent  flow  through  the  pore  (Jvp)  is  obtained  by  integrating  each 
velocity  term  over  the  appropriate  pore  area: 


Jvp  =  Jvp2itrdr  +  J  vn2 itrdr 


(65) 


Solving  these  integrals  and  grouping  terms  yields  an  expression  for  solvent 
flux  in  terms  of  pressure  gradients,  and  pore  and  solvent  radii: 


dx 


(1 


2dU 

dx 


(66) 


where  X=  rs/rp.  Equation  66  may  be  expressed  using  the  coefficient  1^ 
and  a  as  follows: 


Jy  =  kyW 


dP  dU 

—  -  a  — 

dx  dx 


where: 


k 


V 


8tiw 


(67) 


(68) 


o  ={1-(1-X)2)2 


(69) 


The  volumetric  flux  equation,  extended  to  a  solution  of  multiple,  non- 
interacting  solutes,  is  presented  in  Box  5. 


Equations  describing  solute  flux  through  capillary  pores  use  the 
same  geometric  model  as  that  used  to  describe  solvent  flux  (Figure  1).  For 
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simplicity,  total  solute  flux  is  divided  into  pure-diffusion  and  pure- 
convection  components.  To  determine  the  rate  of  pure  solute  diffusion 
(assuming  solvent  flux  is  zero),  a  small  cross  section  of  the  pore,  of 
thickness  dx,  is  considered6  (Figure  2).  Because  the  centers  of  the  solute 
molecules  are  restricted  from  the  area  close  to  the  pore  walls,  the  effective 
area  for  diffusion  (A^,)  includes  the  central  area  of  the  pore  out  to  a 
radius  of  rp  -  rs.  The  cross-section  is  considered  to  be  bounded  upstream 
and  downstream  by  semipermeable  membranes.  A  difference  in  solute 
concentration  across  the  section  leads  to  a  difference  in  osmotic  pressure 
(dll).  If  unopposed,  this  osmotic  pressure  gradient  causes  water  to  flow' 
from  one  side  of  the  section  to  the  other.  In  a  closed  system,  this  osmotic 
pressure  gradient  is  associated  with  movement  of  the  cross  section  in  the 
opposite  direction.  To  prevent  movement  of  the  segment  and  to  maintain 
the  osmotic  pressure  gradient,  a  force  must  be  exerted  on  the  segment 
contents  to  prevent  osmotic  flow.  Thus,  each  solute  molecule  w-ithin  the 
section  is  considered  to  be  acted  upon  by  a  pressure  force  (Fn  )  of 
magnitude: 


dU 

NACsdx 


(70) 


where  the  solute  concentration  within  the  section  is  Cs,  and  NA  is 
Avagadro’s  Number.  In  a  system  with  constant  velocity,  this  force  must  be 
balanced  by  the  drag  on  the  solute  molecule  caused  by  solute-solvent 
interactions  and  interaction  with  the  pore  wall.  Faxen14  determined  that 
the  appropriate  drag  force,  which  accounts  for  the  presence  of  the  pore 
wall  is: 


-  6m\  rsvs 


f. 


(71) 


where  -6xrjrsvs  is  the  Stoke’s  drag  force  for  a  sphere  in  an  unbounded  fluid, 
and  fsp  is  a  correction  coefficient  for  the  presence  of  the  pore  walls: 

f  =  1  -  2.104433A  +  2.08877 A3  -  0.94813X5  -  1.372X6  ♦  3.87A*  -  4.19X10 

(72) 

At  constant  velocity,  Fn  equals  FD.  Setting  the  right  side  of  Equation  70 
equal  to  the  right  side  of  Equation  71  and  solving  for  vs  yields  an 
expression  for  solute  velocity  when  volumetric  flow  is  zero  (i.e.,  pure 
diffusion): 
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VW  " 


[NACt6n^j 


dn 

dx 


(73) 


Total  solute  flux  due  to  diffusion  for  the  single  pore  (Jscjllf)  is  obtained  by 
multiplying  vs<Jjff  by  the  number  of  moles  of  solute  at  the  pore  cross-section 
(Cs  x(rp-rs)2).  After  simplification,  Jsdi(r  is  expressed  by: 


•W  * 


-fs 


{NA6xr\rs) 


dU 

dx 


(74) 


To  determine  the  rate  of  pure  convective  solute  transport  (assuming 
concentration  gradients  equal  zero)  in  a  single  pore,  we  assume  that  the 
center  of  a  neutrally  buoyant  sphere,  suspended  in  a  Poiseuillian  flow, 
translates  along  a  streamline  of  the  undisturbed  flow15.  The  velocity  of 
such  a  solute  molecule,  located  at  a  distance  r  from  the  pore  axis,  is  the 
same  as  the  Poiseuille  fluid  velocity  (vf)  at  that  point: 


VmtvW  “  v/(r)  =  2 


(75) 


where  vf  is  the  average  fluid  velocity  within  the  pore  cross-section. 
Convective  solute  flux  for  the  pore  is  obtained  by  integrating  solute  velocity 
over  the  pore  area  bounded  by  r=0  and  r=rp-rs,  and  multiplying  by  solute 
concentration: 


Jr~T> 

J,,co*v  =  Cs  J  V2,tr<fr 
0 


(76) 


Substituting  the  right  side  of  Equation  75  for  vs  conv,  solving  the  integral  and 
simplifying  terms  yields  an  expression  for  the  convective  solute  transport  in 
a  single  pore: 


-  2 


s.conv 


Csvf*r;(l-Xf[2-(1-Xf] 


(77) 


where,  again,  X  is  rs/rp. 
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Total  solute  flux  (convection  plus  diffusion)  for  a  single  pore  is 
obtained  by  adding  Equations  74  and  77,  substituting  for  the  term  it(rp- 
rs),  and  substituting  Jv  for  the  term  vf7rrp2: 


Js,p 

dn 

dx 

*  Cs4(1-1)2[2-(1-A)j) 


(78) 


Equation  78  may  be  rewritten  in  a  form  similar  to  the  non-ideal  solute 
transport  equation  (Equation  52)  by  regrouping  terms: 


-  CJW 


dU 

dx 


+  %C  J 


where: 


NA6nr\rsw) 


(79) 


(80) 


X  =  (1  “  *)2[2  -  (1  -  X)2]  <81> 

The  solute  transport  equation,  for  a  solution  of  multiple  solutes,  is 
presented  in  Box  5, 


MODEL  COMPARISONS 


The  choice  of  model  equations  for  transcapillary  flux  depends  upon 
the  application.  Although  the  Kedem-Katchalsky  model  has  been  criticized 
as  being  too  simplistic,  this  criticism  may  not  apply  to  studies  of  whole- 
body  net  transcapillaiy  flux.  To  assess  the  effects  of  different  model 
assumptions  on  predictions  of  plasma  volume  response  to  hemorrhage  and 
fluid  resuscitation,  a  simple  two  compartment  model  was  devised. 

Constants  and  initial  conditions  for  this  model  are  provided  in  Table  1. 
Model  compartments,  representing  the  vascular  and  interstitial  fluid  spaces, 
were  considered  compliant  and  well-mixed.  Hydrostatic  pressures  were 
estimated  from  compartment  volumes  as  follows: 
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?c  =  0-91  (VB  -  V^)  +  PCa 

(82) 

Pj  -  0.01042(F)-  VjJ  +  Pjo 

(83) 

where  VB  is  the  total  blood  volume,  V,  is  interstitial  volume,  and  the 
subscript  0  denotes  initial  values.  For  simplicity,  it  was  assumed  that 
plasma  and  interstitial  fluids  were  made  up  of  water,  NaCl,  albumin,  and 
globulin.  At  body  temperature,  the  osmotic  pressure  for  NaCl  (subscript  s) 
is  approximately: 

H,  =  (31. 877)  ( 19330)  Cs  (84) 

where  Cs  is  the  concentration  of  NaCl  (g/ml),  31.877  converts  grams  of 
NaCl  to  milliosmoles,  and  19.33X103  is  the  osmotic  pressure  (mm  Hg) 
exerted  by  1  wosm/ml  of  solute16.  Cubic  equations  presented  by  Landis 
and  Pappenheimer17  are  used  to  estimate  partial  osmotic  pressures  for 
albumin  (subscript  a)  and  globulin  (subscript  g): 

Ea  =  280 Ca  +  18Q0C2  +  12000C3  (85> 

=  210(Cfl  +  Cg)  *  1600 (Ca  +  Cg)2  +  9000 (C0  +  Cgf  -  na  (86) 

where  Ca  is  albumin  concentration  (g/ml)  and  Cg  is  globulin  concentration 
(g/ml). 


Fluid,  NaCl,  and  albumin  exchange  between  the  blood  and 
interstitial  compartments  was  assumed  to  occur  across  the  capillary  wall,  as 
well  as  by  lymph  return  (transport  of  globulin  was  neglected).  Parameter 
values  for  the  flux  equations  were  obtained  by  requiring  that  at  normal 
steady-state,  total  volumetric  flux  is  equivalent  to  lymph  return  (0.02  ml* 
kg'1  ^min'1),  and  that  total  solute  flux  is  equivalent  to  convective  solute 
transport  in  lymph  (interstitial  concentration  multiplied  by  normal  lymph 
flow  rate).  In  addition,  for  models  based  on  irreversible  thermodynamics, 
solute  reflection  coefficients  were  assumed  to  equal  the  values  provided  in 
Table  2.  Initial  parameter  values  for  the  pore  theory  model  were  based  on 
reported  values  for  solute18  and  pore19  radii  identified  in  Table  3. 

Because  pore  theory  parameter  values  were  not  comparable  to 
parameter  values  for  the  models  based  on  irreversible  thermodynamics,  it 
was  necessary  to  run  two  sets  of  comparisons.  In  the  first  set,  we 
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compared  predictions  generated  by  the  three  models  based  on  irreversible 
thermodynamics.  There  was  no  appreciable  difference  between  predictions 
generated  by  these  three  models  in  response  to  a  simulated  instantaneous 
a)  25%  loss  of  blood  volume,  b)  50%  increase  in  plasma  NaCl 
concentration,  c)  50%  increase  in  plasma  albumin  concentration,  or  d)  25% 
loss  of  blood  volume  with  a  50%  increase  in  plasma  NaCl  concentration 
(Figure  3).  Therefore,  it  was  concluded  that  considerations  of  membrane 
continuity  and  solution  non-ideality  are  unimportant  in  predictions  of 
whole-body  plasma  volume  response  to  changes  in  blood  volume  or  plasma 
osmolarity. 

The  second  set  of  comparisons  included  a  discrete,  ideal  form  of  a 
two-pore  hydrodynamic  model  (Box  6)  and  a  two-pore  parameter-averaged 
model,  identical  in  form  to  the  Kedem-Katchalsky  model  (Box  2). 

Parameter  values  for  the  two-pore  model  are  provided  in  Table  3. 

Averaged  parameter  values  are  provided  in  Table  4.  The  two  models 
generated  nearly  identical  predictions  of  plasma  volume  for  simulation 
conditions  a)  through  d),  above  (Figure  4).  These  results  indicate  that  the 
effects  of  capillary  wall  nonhomogeneity  may  be  accounted  for  using  a 
simple,  parameter-averaged,  phenomenologic  model. 

The  response  of  the  two-pore  model  to  simulated  changes  in  plasma 
solute  concentration  or  total  blood  volume  (Figure  4)  is  much  slower  than 
that  predicted  by  models  based  on  irreversible  thermodynamics  (Figure  3) 
or  that  observed  under  laboratory  conditions20.  Modification  of  pore  sizes 
(rsm<2.3  A,  rlg=60  A)  enables  the  two-pore  model  to  achieve  a  slightly 
faster  response.  However,  this  response  is  still  much  slower  than  that 
predicted  using  models  based  on  irreversible  thermodynamics.  It  is 
possible  that  extending  the  pore  theory  model  to  three  or  more  different 
pore  sizes  will  enable  the  model  to  achieve  more  reasonable  plasma 
volume  estimates.  The  models  based  on  irreversible  thermodynamics 
already  predict  reasonable  plasma  volume  responses,  however,  so  the 
complexity  of  a  3-or-more-pore  model  does  not  appear  to  be  warranted. 
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SUMMARY 


1  ae  Kedem-Katchalsky  equations3  make  three  assumptions  which 
might  restrict  their  use  for  describing  transcapiliary  flux:  1)  that  the 
membrane  is  sufficiently  thin  to  express  concentration  gradients  in  discrete 
form,  2)  that  the  solutions  on  either  side  of  the  membrane  are  ideal,  and 
3)  that  the  capillary  wall  is  homogeneous.  Several  investigators  have 
attempted  to  remedy  these  limitations  by  reworking  the  Kedem-Katchalsky 
equations  for  the  more  general  case  of  a  relatively  thick  membrane  (Patlak 
model4),  non-ideal  solutions  (Katz  model5),  and  capillary  wall 
nonhomogeneity  (pore  theory  models6).  The  results  of  comparisons  in  the 
previous  section  show,  however,  that  plasma  volume  predictions,  generated 
by  models  incorporating  different  transcapiliary  flux  equations,  were 
indistinguishable.  Thus,  the  Kedem-Katchalsky  equations,  when 
incorporated  into  a  simple  two-compartment  model,  are  as  effective  as 
more  complex  models  for  describing  the  plasma  volume  response  to 
hemorrhage  and  resuscitative  fluid  administration. 
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Table  1.  Initial  conditions  and  parameter  values  for  2-compartment 
(vascular  and  interstitial)  model. 


blood  compartment 

Blood  Volume  (BV) 

70  ml/kg 

Plasma  Volume  (PV) 

38.5  ml/kg 

Plasma  NaCl  Concentration  (Ct  ? ) 

0.0093  g/ml 

Plasma  Albumin  Concentration  (CaP ) 

0.045  g/ml 

Plasma  Globulin  Concentration  (Cg  P ) 

0.025  g/ml 

Capillary  Hydrostatic  Pressure  (Pc ) 

20  mm  Hg 

interstitial-tissue  compartment 

Interstitial  Volume  (IV) 

140  ml/kg 

Tissue  Volume  (TV) 

350  ml/kg 

Interstitial  NaCl  Concentration  (Cs , ) 

0.0093  g/ml 

Interstitial  Albumin  Concentration  (Ca j ) 

0.016  g/ml 

Interstitial  Globulin  Concentration  (C&! ) 

0.0  g/ml 

Interstitial  Hydrostatic  Pressure  (Pj ) 

-0.223  mm  Hg 
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Table  2.  Parameter  values  for  models  based  on  irreversible 
thermodynamics. 


assumed  values 


K 


0.109 

0.94 

1.0 

35.71  ml*  kg'1 'min'1 


/ 


100x1  O'8  m 


calculated  values 

ka 


0.0612  ml* kg  ’ •min1* mm  Hg'1 
0.0103724  ml ‘kg'1*  min'1 
2.238  x  105  g-kg^min’-mm  Hg 
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Table  3.  Pore  theory  model  parameter  values  for  a  system  with  two  pore 
sizes,  large  (subscript  lg)  and  small  (subscript  sm). 


40X10-8  m 
200x1 0"8  m 
lOxlO^m 
36X10-8  m 
56x10"®  m 


2.252  x  10‘9  ml -kg'1- min'1 
9.970  x  10'7  ml -kg1*  min'1 
5.602  x  1012  ml -kg'1  •min'1 
1 .409  x  10'10  ml  •  kg'1  •  min'1 
0.0  ml ‘kg'1 'min'1 
4.0 183x1  O'12  ml -kg'1 -min'1 
0.01055 
0.0005229 
0.9401 
0.1073 
1.0 

0.9894 

0.9995 

0.0598 

0.8927 

11270160 
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Table  4.  Flow-averaged  parameter  values  for  the  one-pore  model 


K 

0.0264855  ml* kg*5 -min"1 

K 

6.32948x10s  ml* kg*1* min*5 

K 

4.464988x1  O'9  ml* kg*5 -min*1 

0.0103229 

a„ 

0.9053186 

1.0 

.  8 _ 
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Box  1.  Landis-Starling  Equation 


fixation 


'  kl<pc  -  f,)  -  (He  -  n/)j 


Membrane  is: 

1)  Freely  permeable  to  small  solutes 

2)  Non-permeable  to  proteins. 

II  represents  protein  oncotic  pressure 


Box  2.  Kedem-Katchalsky  Equations  extended  to  a  system  of  multiple, 
non-interacting  solutes. 


assumptions 


j,  -  -  £  °.An»j 


J,  -  k*c,  *  a-o.w.c. 


Constant  T  and  P 
Negligible  solute  volume  fractions 
Thin  membrane 
Ideal  solutions 

Constant  partial  molar  volumes 


J  Aln(C,) 
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Box  3.  Patlak  Equations  extended  to  a  system  of  multiple,  non-interacting 
solutes. 


equations 

Jv  =  *V[AP  -  £  c,ADj 

dC 

J,  -  +  ci  -o,)Jvc, 


assumptions 

Constant  T  and  P 

Negligible  solute  volume  fractions 

Ideal  solutions 

Constant  partial  molar  volumes 

Homogeneous  membrane 

Lp,  ws,  and  as  are  independent  of  Cs 


Box  4.  Continuous  model  with  considerations  for  non-ideality  extended  to 
a  system  of  multiple,  non-interacting  solutes. 


equations 

assumptions 

Jv  =  kv [AP  -  JXAn,] 

Constant  T  and  P 

Constant  partial  molar  volumes 
Homogeneous  membrane 

J  =  -  0)  w -  +  (l-a  )J  C 

3  3  V  3'  V  J 
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Box  5.  Basic  pore  theory  model  extended  to  a  system  of  multiple,  non¬ 
interacting  solutes. 


assumptions 


Non-interacting  solute  molecules 
Navier-Stokes  flow 
Cylindrical  pore 

Solutes  do  not  affect  solvent  flow 


equations 


-  kaap  T0./11,] 


uu 

J  =  -  (o  w  — £  +  y  C  J 

StP  StP  1.  A‘X.DX's*'\ 


<3 

fa 


k  = 

v'p  8i!w 


2\2 


- 11  -O-V2) 


„ 


/«  1 


HvMl-V'P-d-A,/] 
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Box  6.  Discrete,  ideal  pore  theory  model. 


Non-interacting  solute  molecules 
Navier-Stokes  flow 
Cylindrical  pore 

Solutes  do  not  affect  solvent  flow 
Ideal  solutions 
Thin  membrane 
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on  irreversible  thermodynamics. 
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